Evolution of coupled classical fields 
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We study the evolution of the coupled scalar and fermion fields within the classical field theory. 
We examine the case of N coupled fields in 1 + 3 dimensional space. The general expressions for 
the fields distributions are obtained. The particular case of two fields in 1 + 1 dimensional space 
is carefully studied. We obtain the expressions for the averaged fields intensities and show that 
in the relativistic limit they are similar to the usual transition probabilities formulae of neutrino 
oscillations. 
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The particle mixing plays an important role in elemen- 
tary particle physics. According to the experimental data 
the particle mixing exists in both quark and lepton sec- 
tors of the standard model. The idea of mixing among the 
two quark flavors was put forward in Ref. [l| to explain 
the baryons decays. The mixing in the leptonic sector 
of the standard model was proposed in Ref. . In that 
parer the neutrino mixing and oscillations were studied 
on the analogy of the known at that time K° <-> K° oscil- 
lations. Then in Ref. this approach was generalized on 
the mixing between three neutrino flavors. Recently we 
obtained the strong evidence in favor of neutrino oscilla- 
tions and, therefore, mixing (see, for instance, Ref. 4]). 
For example, neutrino oscillations are likely to be the 
most plausible explanation of the solar and atmospheric 
neutrino problems. 

In Ref. the neutrino oscillations were examined 
within the quantum mechanical approach. Schrodinger 
like differential equation for the description of the two- 
level neutrino system was proposed. On the basis of this 
equation one can derive the famous transition probability 
formula, 



P(t) =sin 2 (26» vac )sin 



Am 2 



(1) 



where # V ac is the vacuum mixing angle, Am 2 is the mass 
squared difference and E is the energy of the system. Up 
to now Eq. Q is of use in numerous phenomenological 
studies of neutrino oscillations. However more profound 
analysis of particle mixing and oscillations is necessary. 
The approach to the oscillations phenomenon based on 
the field theory methods should be elaborated. 

In the last decade a great deal of studies on the 
field theoretical substantiation of Eq. (JJJ were carried 
out. First of all it is necessary to mention works by 
M. Blasone and G. Vitiello and their collaborators (see 
Refs. HB0H). In these papers the authors made the 
comprehensive analysis of the fermion and boson mixing 
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transformations using the methods of quantum field the- 
ory. It was demonstrated that the vacuum structure of 
mass eigenstates is not equivalent to one of flavor eigen- 
states. The quantum mechanical formula for the tran- 
sition probability was also reproduced. Moreover some 
corrections, which result from more careful quantum field 
theory analysis, were obtained. In Ref. [9( the group the- 
oretical aspects of neutrino oscillations were discussed. 
Analogous approach to description of the neutrino oscil- 
lations was developed in Refs. 03[ll|. 

Rather appreciable contributions to the investigation 
of the flavor neutrino oscillations were made in Refs fl^ . 
13] . In those papers the wave packages treatment of neu- 
trino oscillations was developed as well as the discussion 
of the neutrino oscillations phase is presented. The neu- 
trino oscillations phase was also studied in Ref. 0]. A 
very interesting approach to the description of the neu- 
trino flavor oscillations was proposed in Ref. |15| . In that 
paper the covariant path amplitudes method was applied 
for the analysis of the neutrino oscillations phase. 

Recently we elaborated the quasi-classical approach for 
the description of spin (see Re fe. fltl ITtIIT^I ) and flavor 
neutrino oscillations (see Ref. [l!j ). It was shown that 
neutrino oscillations in moving and polarized matter un- 
der the influence of arbitrary electromagnetic fields were 
described by the generalized Lorentz invariant quasi- 
classical Bargmann-Michel-Telegdi equation. It is inter- 
esting to note that the equation describing the precession 
of the neutrino three dimensional spin vector (neutrino 
spin in particle's rest frame) is the usual Bloch equation. 
We demonstrated that neutrino spin rotates around a 
certain direction determined by the velocities and polar- 
izations of background fermions as well as the electro- 
magnetic field strength. 

The method involving the Bloch equation for the treat- 
ment of neutrino flavor oscillations was proposed in 
Ref. . If one considers the evolution of two neutri- 
nos system (e.g., v e and v^), it is possible to introduce 
the "polarization" vector P = Tr(<rp), where er are the 
Pauli matrices and p is the 2x2 density matrix. If 
neutrinos propagate in vacuum, the vector P was shown 
to precess without loss of length according to the Bloch 
equation. The appearance of classical effects in various 
quantum systems (including the analysis of a two level 
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system with help of the Bloch equation) was discussed in 
Ref. [2l| . Thus basing on the similarity in the description 
of neutrino spin and flavor oscillations we suppose that 
classical theory methods could have been applied for the 
treatment of the flavor oscillations. However this suppo- 
sition should be substantiated by the direct calculations 
that show the classical theory yields at least the same 
results as the quantum one. 

In this paper we study the evolution of the coupled 
scalar as well as fermion fields within the context of classi- 
cal field theory. The main goal of our article is to demon- 
strate that neutrino oscillations can be described within 
the classical approach. The classical approach was also 
adopted since we should not be puzzled by a problem: 
must we rely on flavor or mass eigenstates in our treat- 
ment of neutrino oscillations? The intensive discussion 
about this topic takes place in Refs. 0,H3|. The case 
of N coupled fields in f + 3 dimensional space is exam- 
ined. We solve the Cauchy problem for this system, i.e. 
for the given initial conditions we find the fields distribu- 
tions for any time point. In order to analyze the obtained 
expressions we study the particular case of two fields in 
1 + 1 dimensional space. For the specific initial conditions 
the expressions for the averaged fields intensities are ob- 
tained. We also show that in the relativistic limit they 
are similar to the usual transition probabilities formu- 
lae of neutrino oscillations in vacuum. It is interesting 
to mention that the expressions for the averaged fields 
intensities are identical for both bosons and fermions. 

First let us discuss the case of N arbitrary coupled 
scalar fields. For simplicity we suppose that these fields 
are the real ones. The Lagrangian for this system is ex- 
pressed in the following form 

JV JV 

= ^2^o(v>k) + 9ikfifk, (2) 

k—1 — 1 

i^k 

where gik are the coupling constants, ip — {(fx, . . . , </?jv), 
and 

£o(Vfe) = \d^kd^ k - (3) 

is the Lagrangian for the field ipk(r,t) at the absence of 
the additional coupling, is the mass corresponding to 
this field. It is necessary to note that the second term 
in Eq. © is assumed to be an interaction between fields 
ipk- 

In order to describe the evolution of the system ©-10 
we should set the Cauchy problem for this system. For 
the initial conditions, 

Pi(r,0) = /i(r), ft(r,0)= ffi (r), (4) 

where fi(r) and <?i(r) are the given functions, one should 
find the fields distributions ipk{v,t) for any time point. 
It is always possible to diagonalize the Lagrangian @ 



with help of the transformation, 

JV 

<Pi(r,t) = Y M ikUk(r,t). 
fc=i 

Thus the Lagrangian expressed in terms of the fields 
ttfc(r, t) takes the form 

JV 
k=l 

where Co(uk) is the Lagrangian for the field itfc(r, i), 

A)(wfc) = ^df.Ukd^Uk - -yu k , 

and mfe are the corresponding masses. It should be noted 
that these masses differ from the masses of the fields ipk- 
The fields Uk(v,t) are usually called mass eigenstates in 
contrast to <pk(v,i). 

One can write the differential equations for the fields 
tife(r, t). It is the system of the usual homogeneous Klein- 
Gordon equations. Their solutions have the form 

Uk(r,t) = J ^[a+(p)e-^+a-(p)e^]e^ (5) 

where 8k = \/p 2 + rri l k , and a^(p) are the Fourier coef- 
ficients. Note that ajf (p) are the c-numbers. 

To solve the Cauchy problem we introduce the func- 
tions 

JV 

F fc (r)=^(M- 1 ) H / i (r), 

i=l 

JV 

G fc (r)=^(M- 1 ) fei5l (r). 

i=l 

These functions are the initial conditions for the fields 
life. Then one should pick out the coefficients a^(p) so 
that to satisfy the initial conditions ©. From Eqs. Q 
and JSJ we obtain 

4(P) = ^(P)±^)> 

where -Ffc(p) and Gfe(p) are the Fourier transforms of 
the functions -Ffc(r) and Gfe(r) respectively. Finally we 
receive the fields distributions ipj(v, t) in the explicit form 

JV 

<Pi(T,t) = Y,M ik {M- l ) H 

ik=l 

x J dV[4(r-r',t)/i(r') 

+ £> fe (r-r',t)< 7i (r')], (6) 

where 



3 



is the Pauli-Jordan function. It is interesting to list some 
of the properties of the Pauli-Jordan function 

D k (r,0)=0, D k (r,0)=d 3 (v), D fc (r,0)=0. 

It is worth noticing that the initial conditions in Eq. (@J 
are consistent with these properties of the Pauli-Jordan 
function. We also mention that the Pauli-Jordan func- 
tion can be expressed in the explicit form (see, e.g., 
Ref. H), 



D k (r,t) = ±e(t)S(s 2 ) 

where s 2 — t 2 — r 2 and 

'l, t>0, 
-1, <<0. 
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e(i)0(*Vi ("»**)» (8) 



e(t) = 



[s) = 



1, s > 0, 
0, s < 0. 



are the step functions. Thus Eqs. JSJ and (JJJ) represent 
the exact solution of the Cauchy problem for arbitrary 
functions /j(r) and <?j(r). 

The integrals calculation in Eq. @, however, are 
rather awkward in general 1 + 3 dimensional space. Thus 
let us, for simplicity, consider the space with 1 + 1 di- 
mensions. Instead of Eq. J7J) we have 



D k (x,t) 



+00 



d P jpx sm£ k t 



2tt 



Now £ k = \/p 2 + m 2 . One can also obtain the Pauli- 
Jordan function in the explicit form [see Eq. (JSJj in 1 + 1 
dimensional space, 



D k (x,t) = -9(s 2 )J {m k s), (9) 
D k (x, t) = tS(s 2 ) - ^9(s 2 ) Ji(m feS ). (10) 

Here s 2 — t 2 — x 2 . 

We suppose that q k (x) = and f k (x) ^ 0. Then we 
receive (see also Ref. |24j) 

JV 



y j {x,t) = Y,M jk {M- l ) kt 
=i 

{\[h(x-t) + h{x + t)\ 

2 Jx-t s > 



ik=l 
X 



(11) 



It should be noted that if functions fi(x) ^ in a 
bounded region and f%(x) — > 0, when a; — > ±oo then 
the second term in Eq. (1111) is the vanishing one and 

^(M)^o[/i(a-*) + /i(s + *)]- 



interesting physical implication. If a single particle ap- 
pears far from a detector, then its field distribution is 
localized in space. When a particle begins propagating 
towards a detector its field distribution approaches to the 
initial conditions. Thus the effect of various non-trivial 
phenomena (like conversion, or oscillations, from one field 
type to another) will be vanishing. 

Now let us choose the initial conditions. We suppose 
fi(x) = and 



/ 2 (o0=2lsin(|z), a: 



where L is the "volume" of the space. Note that 21 is just 
the normalization factor. In this case we can calculate 
the integral in Eq. I|ll|l explicitly 



/J dy sin(| 



lu \ Ji(ms) 

y) — : — 



™ (~x) J 1/2 J 1/2 , (12) 



where 



In computation of the integral in Eq. H12(l we used the 
expression 



' cos(6\/a 2 — x 2 ) T . , 

ax -j^=^= — J v {cx) 

V« 2 — x 2 



\Ju/2 [I (VW + #-b)] J v/2 [| + 6) 



One half order Bessel function can be expressed in terms 
of the elementary function. Namely, 



(13) 



J\/n{z) = \ — s i nz - 

" TTZ 



Let us consider the case when u> has great values com- 
pared to the masses m^: cj S> mi,2- This situation 
corresponds to the high energy approximation or rela- 
tivistic "particles". Then, the parameters (3\^ take the 
form 



01 



m 



CO, 



(14) 



The field distribution tpi (x, 0) is equal to zero. Thus we 
can describe the dynamics of this field for the subsequent 
points of time. If one studies the evolution of two fields, 
ipi^ix^t), it is possible to parameterize the matrix Mj k 
with help of one angle 



This property of the Klein-Gordon equation was also 
mentioned in Ref. 24]. The described feature has one 



M jk = 



cos 8 sin 9 



sin v cos ( 



(15) 
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Using Eqs. l|T2 )l -l|T5 fl we can rewrite Eq. in the fol- 
lowing way 



(x,t) =2asin20sin (~x) sin (^tj 



x sm 



x cos 



Auj 



{ml - m\) 



— (m? + mf) 



(16) 



Now let us discuss the field measurement process. In 
case of rapidly varying fields (uj ^> mi :2 ), a detector 
registers not the field strength, but the intensity of the 
field which is proportional to the field strength squared, 
I ~ (p 2 (x,t). Moreover a detector has limited sensitivity, 
i.e. it cannot register arbitrary field variation in time and 
in space. Thus we should average the intensity over the 
characteristic time and space scales. These scales should 
be greater than typical time and space scales of the field 
in question, i.e. 1/uj. 

To calculate the mean value of the intensity one should 
take into account the expressions 



sm —x 
2 



= ( sm' 



and 



t 

Auj 
t 

Aw 



(ml 



(ml + m 2 ) 



t 

AlO 

t 

Auj 



(ml - ml) 
(ml + m\) 



since (mf ± m|)/o; -C uj. Then we should introduce the 
normalized intensity of the field y>i(x, t) according to the 
formula 



P(t) 



2t 2 ' 



Finally we obtain the expression for P(t) in the following 
form, 



P(t) = sin 2 29 sin 2 



AlO 



. 2 / ml + m\ 
x 4 1 - sm I — -t 

Auj 



(17) 



where we introduced the common notation Am 2 = ml — 
m\. 

It is necessary to identify the uj parameter. We cannot 
directly equate it to the particle energy, E = hu, since 
we are using the classical approach here. Moreover, the 
chosen "wave function", /a(^) ~ sm(ujx/2), does not cor- 
respond to a definite momentum and thus to a definite 
energy. However we can calculate the averaged energy 
density of the system, 



(pi 




+ mlfZ(x) 



(18) 



Here we suppose that fi(x) = 0. Using Eq. l(TH|) in rela- 
tivistic limit (u> 3> mi^) we obtain that 



(Pe) 



and we can identify uj with the energy of the system. 
Thus the first term in Eq. I|17(l is similar to the well 
known formula for the transition probability in the two 
neutrino system. It is interesting to mention that the sec- 
ond term in Eq. (|17fl contains the harmonic oscillations 
with the frequency (m 2 



raj! [Auj. Analogous additional 
|E IE Hi an d was treated as 



term was obtained in Refs. 
the quantum field theory correction to the Eq. . How- 
ever our approach demonstrates that this term appears 
when one uses classical field theory. It results from the 
accurate account of the Lorentz invariance. 

Now let us discuss the case of the coupled fermion 
fields. The Lagrangian for this system is expressed in 
the following way 



N 



( 



C(u) = Y,C Q (v k ) 



fc=i 



N 



^ 9%kV%Vk + h.c. , 



where v = (y\, . . . , fjsr), and 

£o(vk) = Vkiil^dn " t"fe)i/fc. 

We again should set the Cauchy problem for the system of 
differential equations in question. However here one has 
to impose only one initial condition since Dirac equation 
is the first-order differential equation, 



^(r,0)=6fc(r). 



(19) 



Analogously to the case of the scalar fields we can intro- 
duce the mass eigenstates, 



N 



vfat) = > M ik il> k (r,t). 



fe=i 



The Lagrangian expressed in terms of the mass eigen- 
states has the form, 



JV 



£(TA) = ^£ (V'fe) 1 



fc=i 



where 



is the Lagrangian for the field V'fe( r i t). Note that masses 
mfc differ from the masses of the fields Vk ■ 

The solution of the Dirac equations for the mass eigen- 
states fields can be expressed in the following way, 



ipk 



fe(r ' t)= /(2^)^[ as(p)Ms(p) 



-iSkt 



+ b s (-p)v s (- P y £kt ]e* pr . 



(20) 
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Here u s (p) and v s (p) are the basis spinors, a s (p) and 
6 s (p) are the indeterminate functions. 

Now we should find the values of the a s (p) and b s (p) 
c-number functions to satisfy the initial condition given 
in Eq. I|19l) . The calculations are analogous to the previ- 
ously discussed case of the coupled scalar fields. Thus we 
arrive to the solution of the Dirac equations which are 
valid for arbitrary functions £j(r), 



A' 



ik=l 



x / d 3 r'S fc (r'-r,t)(-i7 )6(r'), (21) 



where 



S k (r,t) = {ij^+m^Dk^t), x» = (t,r), 

is the Pauli- Jordan function for a fermion field (see, e.g., 
Ref. In deriving of Eq. we used the orthonor- 

mality conditions 

<4(p)«r(p) = "I(pH(p) = $sr, 

*4(pR(-p) = vl(p)u r (-p) = 0, 
and the formulae for the summation over the spin indexes 

I>.(k)«t{k) = ^ 7 °, 

x;«.(k)«t(k)=^^ 7 °. 

It is interesting to mention that the function S k (r, t) has 
the following property, 

S k (r,0) = i 1 °5 3 (r). 

Thus the solution given in Eq. Q21|l is consistent with 
the initial condition (|19|) . Equation (|21[1 can be rewrit- 
ten in the non-covariant form which is, however, more 
convenient for the further analysis, 



N 



zfc=l 
X 



=i 

{-J d 3 r'(aV r )D fe (r-r',i)&(r') 
J dVD fc (r-r',%(r') 
im k (3 f dVD fc (r-r',t)ei(r')}, (22) 



where we use common notations for the gamma matrixes, 
a = 7°7 and (3 = 7 . 

Just for simplicity we again discuss the case of two 
coupled Dirac fields in the space with 1 + 1 dimensions. 
Dirac equation in 1 + 1 dimensional space was carefully 



studied in Refs. |26l 1271 l28j. The gamma matrixes have 
the form, 



7° = 



1 

-1 



7 = 



1 

-1 



(23) 



Now one should set the initial conditions. Let us as- 
sume that £i(x) = and £2(2) is expressed in the follow- 
ing way, 



1 / cos(wa;/2) 

2 \^sin(o;a;/2) 



(24) 



We obtained the expression for the evolution of the 
V\(x, t) which accounts for the exact dependencies on the 
particles masses. However it appeared to be rather awk- 
ward. Nevertheless it can be shown that the the third 
term in Eq. I|22|) is negligible in 1 + 1 dimensional space. 
Indeed, let us consider, for instance, the integral 

I{x,t)=J dyD k (x-y,t)sm(^ y y (25) 

We remind that the Pauli- Jordan function in 1 + 1 dimen- 
sional space is given in Eq. ©. The integral in Eq. I|25|) 
can be calculated explicitly and expressed in the form 



I{x,t) = sin (^xj 



uj \ sm(tyjm 2 k + (w/2) 2 ) 



(26) 



Here we used the known value of the integral 



dx Jo(b\/ a 2 — x 2 ) cos(cx) 



\Jb 2 + c 2 



Thus in the high energy approximation (u> m k ) we 
obtain that I(x, t) — > 0. It is also interesting to note 
that one should carefully follow the order of integration 
and differentiation while using the Pauli- Jordan function. 
Indeed using, for example, Eqs. (|10|) - (|12[l an d (|26|) we can 
see that 



dt 



dxD k (x,t)f(x) ± 



r+00 g 



because Pauli-Jordan function is the singular one. 
Finally we get the expression for the v\{x, t), 



fi(x,t) — sin 28 



x sm 



sin I — tj + cos y—tj 



t 



— TOo) 



cos 



— K + m^) 



/ cos(wa;/2) 



I sin(cja;/2) 

In deriving of Eq. (|27|l we used the fact that 



d_ 

dx 



D k {x,t) 



x d 
"tdi 



D k (x,t). 



(27) 



(28) 
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The measurable quantity of the classical Dirac field is 
the intensity. It is proportional to the \vi(x,t)\ 2 . How- 
ever we again should average the intensity over space and 
time. Thus, using Eq. l|77|l we obtain for (I)(t) the fol- 
lowing expression 

(/>(*) - sin 2 29 sin 2 (^-tj 

^-^(m^')}' (29) 

which coincides with the similar expression derived for 
the scalar field. Note that Eq. (|29|) again contains the 
additional term oscillating with the frequency (mf + 
mf.)/4w. 

The calculations performed in this work demonstrate 
[especially Eq. (|29[) ] that neutrino flavor oscillations can 
be treated in frames of the classical theory. According to 
the classical field theory approach the evolution of flavor 
neutrinos is described in the following way. 

1. Flavor neutrino emission in a reaction. This pro- 
cess can be described by means of the quantum 
approach. However here we should obtain the final 
field distribution rather than the emission proba- 
bility. It is also possible to admit that the mixture 
of the neutrino flavor s ap pears in a process, as it 
was proposed in Ref. |22|. In this case one should 
set other initial conditions in Eqs. (QJ and H19|) . 

2. Neutrino propagation towards a detector. One can 
successfully use the methods elaborated in this pa- 
per for the description of the neutrino conversion 
or oscillations. Basing on the initial fields distri- 
butions obtained in the item ^ we derive the final 
fields distributions that take into particles mixing. 

3. Neutrino interaction with a detector. This process 
again can be described by means of quantum field 



theory. On the basis of the fields distributions ob- 
tained in the item [21 one can calculate the neutrino 
flux measured with a detector. 

Thus we should not directly involve the mass eigenstates 
if we are using classical approach. 

In conclusion we mention that the evolution of the cou- 
pled scalar as well as fermion fields within the classical 
field theory has been studied in this paper. We have 
examined the case of N coupled fields in 1 + 3 dimen- 
sional space. The Cauchy problem has been formulated 
for these systems. We have solved it for arbitrary initial 
conditions. The particular case of two coupled fields in 
1 + 1 dimensional space has been studied. Finally we have 
obtained the expressions for the averaged fields intensi- 
ties. It has been shown that in the relativistic limit these 
expressions were similar to the usual transition probabili- 
ties formulae of neutrino oscillations in vacuum. The dis- 
cussion of the additional terms in transition probabilities 
formulae has been presented. It has been demonstrated 
that the expressions for the averaged fields intensities for 
both bosons and fermions turned out to be identical. We 
have shown by means of the direct calculations that the 
flavor oscillations phenomenon could be described within 
the classical approach. Thus one can conclude that the 
usage of the quantum mechanics is inexpedient because 
classical field theory yields more elegant description of 
the problem in question. 
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